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Abstract 

Within the framework of the minimal non-minimal supersymmetric standard 
model (MNMSSM) with tadpole terms, CP violation effects in the Higgs sector are 
investigated at the one-loop level, where the radiative corrections from the loops of 
the quark and squarks of the third generation are taken into account. Assuming 
that the squark masses are not degenerate, the radiative corrections due to the 
stop and sbottom quarks give rise to CP phases, which trigger the CP violation 
explicitly in the Higgs sector of the MNMSSM. The masses, the branching ratios 
for dominant decay channels, and the total decay widths of the five neutral Higgs 
bosons in the MNMSSM are calculated in the presence of the explicit CP violation. 
The dependence of these quantities on the CP phases is quite recognizable, for given 
parameter values. 
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I. INTRODUCTION 



Any supersymmetric standard models should have broken supersymmetry (SUSY) in 
order to be phenomenologically realistic. An intensively studied method of breaking the 
SUSY is the inclusion of soft breaking terms in the Lagrangian density of the model 
[1]. Phenomenological analyses of the supersymmetric standard model with soft breaking 
terms might become more realistic as well as more interesting if the soft breaking terms 
contain or induce some complex phases, since the presence of the complex phases may 
give rise to CP mixing in the Higgs sector of the model [2] . 

The simplest version of the supersymmetric standard model is the minimal supersym- 
metric standard model (MSSM), which has just two Higgs doublets in its Higgs sector. 
The MSSM at the tree level cannot accommodate explicit nor spontaneous CP violation, 
because any complex phases in the Higgs sector of the MSSM can always be eliminated 
by rotating the Higgs fields. Even at the one-loop level, spontaneous CP violation is 
disfavored by the MSSM because a very light neutral Higgs boson is required, which has 
already been ruled out by experiments [3]. On the other hand, it is known that explicit 
CP violation is viable in the MSSM at the one-loop level, since the radiative corrections 
due to quarks and squarks yield the CP mixings between the neutral Higgs bosons [4-7]. 

A number of nonminimal versions of the supersymmetric standard model have been in- 
troduced in the hterature, by including a Higgs singlet to the Higgs sector of the MSSM [8], 
such as the next-to- minimal supersymmetric standard model (NMSSM) with Z3 symmetry 
[9-16], the general NMSSM with broken Z3 [17], the minimal non-minimal supersymmetric 
standard model (MNMSSM) with tadpole terms [18], and the minimal supersymmetric 
models with an additional U{1) [19,20], to cite a few of them. An advantage of them over 
the MSSM is that the vacuum expectation value of the Higgs singlet may dynamically 
solve the problem of dimensional /i-parameter in the MSSM [21]. 

The NMSSM, which is most popular among those nonminimal versions of the su- 
persymmetric model, has been extensively studied in various aspects of phenomenology, 
including explicit CP violation [13-16]. In the NMSSM, explicit CP violation may take 
place at the tree level, since the NMSSM may have one nontrivial CP phase after re- 
defining the Higgs fields. By assuming the degeneracy of the stop quark masses in the 
Higgs sector of the NMSSM, it is found that large explicit CP violation may be realized 
as the vacuum expectation value of the neutral Higgs singlet in the NMSSM approaches 
to the electroweak scale [13]. We have elsewhere studied the effects of explicit CP viola- 
tion in the NMSSM on the neutral and charged Higgs boson masses at the one-loop level 
by considering radiative corrections due to various particles and their superpartners [15]. 
Also, we have recently calculated the Higgs decays within the context of the NMSSM 
with explicit CP violation [16]. 

The MNMSSM is different from the NMSSM in the sense that explicit CP violation 
is not possible at the tree level in the MNMSSM. This is because the Higgs sector of 
the MNMSSM, unlike the NMSSM, can always absorb the CP-violating complex phase 
by rotating the relevant Higgs fields such that its tree-level Higgs potential cannot have 
any CP phase. Thus, within the context of explicit CP violation in the Higgs sector, the 
MNMSSM is similar to the MSSM with an additional U{1) rather than the NMSSM. At 
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the tree level, the MSSM with an additional U{1) can always eliminate any complex CP 
phase in its Higgs sector by employing the tadpole minimum condition or by rotating the 
three neutral Higgs fields [20]. The CP symmetry of the MNMSSM may be violated at 
the one- loop level in an explicit way as radiative corrections are taken into account. 

In this article, we study the explicit CP violation scenario at the one-loop level in 
the Higgs sector of the MNMSSM. We calculate the radiatively corrected masses of the 
neutral Higgs bosons in the MNMSSM. The dominant decay modes for the neutral Higgs 
bosons into heavy fermion pairs, gluon pairs, and weak boson pairs, are studied. We are 
interested in the dependency of those decays on the CP phases in the explicit CP violation 
scenario at the one-loop level. The branching ratios for the neutral Higgs bosons with 
CP-undefined states are found to vary significantly against the CP phases, arising from 
the squark masses of the third generation. 



II. THE CP- VIOLATING HIGGS POTENTIAL 

The Higgs sector of the MNMSSM consists of two Higgs doublet superfields Hi = H{), 
H2 = (-f^2^, and a Higgs singlet superfield A^. Keeping only the Yukawa couplings for 
the third generation of quarks, the superpotential of MNMSSM may be written as 

W = hteijQ%H'i - h€ijQ%Hl - Ke,jUT%H{ + \eijH\H'iN (1) 

where eij is totally antisymmetric with £12 = —£21 = 1, ^"^^ L are the SU(2) doublet 
quark and lepton superfields of the third generation, respectively, t^, 6^ and f|j are 
the SU(2) singlet top, bottom, and tau superfields respectively, ht, hf,, and hr are the 
Yukawa coupling coefficients of top, bottom, and tau superfields, respectively, and A is a 
dimensionless coupling coefficient. 

The tree-level Higgs potential, V^, of the MNMSSM may be decomposed into D-terms, 
F-terms, the soft terms, and the tadpole terms as 

= Vd + Vf + Vs + Vt , 

where 
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Vt = -ieN + R.c.) , (2) 

with gi and g2 being the U{1) and SU{2) gauge coupling constants, respectively, a be- 
ing the Pauli matrices. Ax being the trilinear soft SUSY breaking parameter with mass 
dimension, niHi, fnH2, and rriN are the soft SUSY breaking masses, and C, is the tadpole 
coefficient. 

Note that the above tree-level Higgs potential would additionally have a global U{1) 
Peccei-Quinn symmetry if there is no tadpole term. The global f/(l) Peccei-Quinn sym- 
metry gives a natural solution to the strong CP problem [22]. However, it leads to the 
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existence of a massless pseudo-Goldstone boson which emerges from the tree-level Higgs 
potential since the determinant of the pseudoscalar Higgs boson mass matrix is zero [23]. 
This is the Weinberg- Wilczek axion. The global U{1) Peccei-Quinn symmetry is eventu- 
ally broken by the quantum effects arising from the triangle anomaly, since the Noether 
current for the global U{1) transformation is anomalous, and the Weinberg- Wilczek axion 
acquires a small mass due to instanton effects. Its mass is inversely proportional to the 
scale of the axion decay constant as 

10^ 

niA = 0.6eV— , 
J A 

where Ja is assumed to be the electroweak scale. By considering the Weinberg- Wilczek 
axion as the candidate for the cold dark matter, cosmology estimates rriA ~ lO""' — 10~^eV 
[24] . This information predicts that the scale of the axion decay constant is bounded as 
/a > 6 X 10^ GeV. 

The Weinberg- Wilczek axion has been excluded by the negative results of experimental 
searches. On the other hand, in such models as the Kim-Shifman-Vainshtein-Zakharov 
model [25] or the Dine-Fischler-Srednicki-Zhitnitskii model [26], the scale of the axion 
decay constant can be very high and thus the axion in these models becomes practically 
invisible. Therefore, it may survive without contradicting the experimental constraints. 

In order to avoid the Weinberg- Wilczek axion in the MNMSSM, we introduce the 
tadpole term in the Higgs potential, as seen in Eq.(2), Thus, the global U{1) Peccei- 
Quinn symmetry in the MNMSSM is explicitly broken. The tadpole coefficient ^ in the 
tadpole term is intrinsically a free parameter. However, it might be assumed that it is 
the same order of the SUSY breaking scale from a phenomenological point of view. 

We would like to show that the Higgs sector of the MNMSSM cannot invoke explicit 
CP violation at the tree level. Note that, if the above tree- level Higgs potential contain 
any complex phases, the possibility would be that Ax, A, or ,^ are complex. Among them, 
without loss of generality. Ax may be assumed to be real, since its phase can always be 
absorbed into the phase of A. Further, A can be made real by adjusting the phases of the 
Higgs doublets, while ^ can also be made real by adjusting the phase of the Higgs singlet. 
Consequently, the Higgs potential of the MNMSSM at the tree level can only have real 
parameters, hence no explicit CP violation. It is also impossible for the Higgs sector of 
the MNMSSM at the tree level to invoke spontaneous CP violation. 

In terms of the physical Higgs fields, the Higgs doublets and the Higgs singlet may 
generally be expressed as 



where Si {i = 1,2,3) are the neutral scalar Higgs fields. Pi {i = 1,2) are the neutral 
pseudoscalar Higgs fields, C"*" is the charged Higgs field, and vi, V2, and x are respectively 




V {V2 + S2 + icospPi)e" 
( X + ^3 + zPa ) e'^ , 



(3) 



N 
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the vacuum expectation values of the neutral Higgs fields with v = + V2 = 174 GeV 
and tan/3 = V2/V1, 9 is the relative phase between Hi and H21 and 5 is the phase of A^. 

Now, the tree-level tadpole minimum conditions upon the tree-level Higgs potential of 
the MNMSSM, which are derived from the first derivatives of the tree-level Higgs potential 
with respect to the two pseudoscalar Higgs fields, may be written as 

= 2\vxA\ sin 61 , 

= Xv'^Ax sin 2/3 sin Oi + 2^^ sin 62 , (4) 

where the two phases 61 and 62 are defined a.s di = 6 + 5 and 62 = S. It is straightforward 
that the tadpole minimum conditions are satisfied only when both 61 and 62 are zero, in 
other words, only if 6 = 6 = 0. Therefore, explicit CP violation is not possible in the 
Higgs sector of the MNMSSM at the tree level. 

In order to accommodate any complex phases so as to induce the CP violation between 
the scalar and pseudoscalar Higgs bosons in the present model, one has to consider higher- 
order corrections. The radiative corrections due to top and stop quarks are known to affect 
significantly the tree-level Higgs sector of supersymmetric models. The contribution of 
the bottom and sbottom quark loops is not negligible for very large tan /3. 

The Higgs potential at the one-loop level may be written as 

V = V'^ + V^ 

where is the one-loop effective Higgs potential including the radiative corrections due 
to quarks and squarks of the third generation, which is explicitly given as [27,18] 

^ =§32^ r^^-2j-T6^r^ A^-2j ' 

where M.q {q = t, h) and M.^^ (i = 1 ,2) are respectively the quark and squark masses 
given as functions of the Higgs fields, and A is the renormalization scale in the modified 
minimal subtraction scheme. 

After the spontaneous breakdown of the electroweak symmetry, the quark masses of 
the third generation are given by = htV2 and = /ifeUi, and the squark masses of the 
third generation are given by 

t2 ~ "^Q + =1= ht\J Afv2 + X^vlx^ + 2\AtViV2X cos (pt , 
^ii h ^ + T 'n^b\/ Al + X^x^ tan^ /3 + 2XAbX tan /3 cos (pb , (6) 

where mg is the soft SUSY breaking mass. At and A}, are the trilinear SUSY breaking 
parameters with mass dimension, and two phases (pt and (pb are given as 

(pt = (pAt + + 6 , 

<Pb = <pAk + + 6 , 

where (pAt and (pA^ are respectively the phases of At and Ab, which are in general assumed 
to be complex. In the expressions for the squark masses. At and Ab are now real parameters 
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because their phases are aheady taken out. Note that the D-terms are not included in 
the squark masses. 

Now, at the one-loop level, the tadpole minimum conditions may be written as 



3hj 



3hl 



= ^Asin^i + ^^Asin</)t/(mj^,mjJ + Y^Asin</)b/(m^^,m^J 

3/^2 

= v4a Aw^ sin 2/3 sin 6*1 + 2^^ sin 6*2 + 777^ At At;^ sin 2/3 sin 0J(m|,m; 



+ 



3hl 



Y^AbXv'^ sin 2/3 sin 0b/ (m?^ , m| ; 



where the scale-dependent function f{ml,ml) is defined as 



(7) 



f{ml,ml) 



ml log — -m^ log^ 



+ 1 



It is quite obvious that = does not satisfy the first tadpole minimum condition 
at the one-loop level. Unlike the tree-level case, Q\ at the one-loop level is not zero but 
dependent on the other parameters. Let us rename Q\ at the one-loop level as 0o hereafter, 
in order to avoid any confusion with the tree-level Q\ that is zero. On the other hand, 
one can easily see that sin ^2 = satisfies the second tadpole minimum condition, by 
substituting the first tadpole minimum condition into the second one. In other words, 
we have ^2 = ^ = and Q\ = 6, renamed as 0o, which is expressed in terms of the other 
parameters. Thus, we are left with (pt and (pi, at the one-loop level. Consequently, the 
Higgs sector of the MNMSSM may eventually have two physical CP phase (pt and (pi, at 
the one-loop level even if there is no complex phase at the tree level. 

The 5x5 symmetric mass matrix, M, for the five neutral Higgs bosons at the one-loop 
level in the MNMSSM is obtained from the second derivatives of V with respect to Si, 
S2, 5*3, Pi and P2- In the basis of {Si, S2, Pi, Ss, P2), where we permute Pi and 6*3 for 
convenience, the matrix elements of M are obtained as 



Mil 
M22 
M33 

M44 

M55 
M12 

Mi3 

Mm 

M23 
M24 
M25 



+ + {mz cos pf + m\ sin^ /?, 
M22 + + {f^z sin + m\ cos^ {3, 



Ml + 



44 



„,2 e3 
— -m^ sm 2p -\ , 

„.2 



X 

3 



M- 
M- 

M- 



^x^ X 
12 + ^12 + (2A^i;^ - ml- m\) sin /3 cos /? , 

rt I A /fb 

V 

—m\ sin^ /3 cos [3 + 2vX^x cos /3 , 

X 



+ Ml - 



Ml + M- 



15 ' 



M23 + M. 



M. 



23 5 



24 



Ml m\ sin /3 cos^ /3 -|- 2vX^x sin (3 , 

X 



Ml + M. 



25 ' 
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M34 = M*4 + M3\, ^ 

M35 = M*5 + Mg'g + ^m^ sin /3 cos /3 , 

M45 = Ml, + Ml, , (9) 

where M*^- (z,j =1-5) are the radiative corrections due to top and stop quarks, M^^ 

= 1-5) are the radiative corrections due to bottom and sbottom quarks, m| = {gl+g2)v'^ /2 

is the squared mass of the neutral gauge boson, and m\ is introduced for convenience as 



., XxAxcos(f)o SmfAtXx cos if. ^ ,-i -i -z 
^ sin /3 cos /3 IGttV sin^ /5 cos/?'' ^ 



j^j/va. L-WD v^j ,2 2 



3ml AbXx COS (pb 2 2 n nn\ 
Ibn'^v^ cos-^ p sm (j "1 "2 



Exphcitly, M* are given as 



Ml 



m: 



M\ 



M. 



M. 



Ml 



SmfA^Al g{m}. mi) ^ SmM^Ar logK/m?) 



Syr^f 2 sin^ /9 (m| — m|)2 4712^2 sin^ /3 (m| — m|) 

3< 1 JK'^i\ 



~l~ o 9 9 • 2/3 1 A I ■ 

Svr^w^ sm p y y 
Sm^A^xM^ sin^ 0t g{mf^, m|) 
~ 87r2?72 sin^ p {m\-m\y 



3m^A^A| ^(m|, m\, 
~ STrHa.n^ (3{m\-m\f ' 
, ^ 3m^AM,^sin^0f m|) 

87r2tan2/5 {m\-m\f ' 
, _ 3m^Ax^Af^A,-^ m|) ^ Sm^AxA,-^ log(m|/m| 



87r2ti2 sin^ /5 (m| — m?J2 87r2t;2 sin^ /5 (m| — m?J 
i ^ _ 3mfA^xMfA,-^sin0t m|) 

«7r2t;2 sin^ /5 (m? - ^2 ' 



^ 3m^A^xA| ^(m|, m|) 3m^A^xcot/3 , 2. 
87r2t;sin/5tan/5(m| -m|)2 Sn'^vsiYip 



.2 ^2 



_ 3m^A2xAAi^sin0t ^?(m|^, 



87r2t)sin/5tan/5 (m? —m^^'^ 



t2 ti' 

?,ml\xAlAl^ sin 0t , m| ) 3m^AxAt sin 0f log(m| /m?^ 



87r2ti2 sin^ /5 (m|— m|)2 8712^2 sin^/3 (m| — m?^ 
3m4AAtAr Ar (7(m?, m?) / 3mnAr log(m?/m?) 



87r2t; sin P tan /3 (m?_^ — m?^ )2 87r2t> sin /3 tan /3 (m?^ — m?^ ^ 



Sm'^XA'^ sin (j)t g{mf^, mf^) 3mfXAtsin(j)t log(m|/m|^ 



87r2t; sin P tan /3 (m2_^ — m?^ )2 87r2t; sin P tan /3 (m?^ — m?^ ' 
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M: 



t _ _ 3m^A^xA Aj-^ sin g{m\, m|) 



Svr^f sin /3 tan /3 (m?^ — m?^ ^ 
^ 3m^A^xA^sin^0, m|) 
3^ 87r2?7sin2/5tan/3 (m| -m|)2 ' 

_ 3m^AMA-,sin0i ^7(m|, m|) 



and Mj^- are given as 



'23 



Svr^f 2 cos^ /3 (m? — m? )2 4-7r2t;2 cos^ /3 (m? — m? 

' ^ 62 01 ' ' ^ 62 01 

. I 61 62 



-I log 

Svr^t^^ cos^ (3 \ ml 

Stt^w^cos^/? (m? — ""^^ )^ ' 
3m4A2xM2 sin^ 0, ^^K^, m?; 



87rVcos4/3 (m?^-m?^)2 



, ^ 3m^A2A2sin^X ^Ki> <^ 

87r2cot2/3 __2 ^2 



I Till — Till 
02 bi ' 



3mtXxA,AjAjQ{ml, ml) ^ 3m^AxA. log(m? /m?) 



62 bl' ^2^ _|_ ^2 °^ ^2^^ bi 

""^^ Svr^t;^ cos2 /? (m? — mr V Stt^v^ cos^ 3 (m? — m? ) ' 

' ^ 02 01 ' ' ^ 62 01 ' 

_ 3m^AxA2A^^sin0fe ^?Kv my 3m^Ax^ sin 0^ log(my m?; 



""^^ 87r2t>2 cos^ /3 (m? — m? )2 STr^t;^ cos^ /3 (m? —m?] 



^ 3m^A^Aj;^Aj,^ ^(my m|) ^ 3m^AA^^ log(m|/my 
""^^ Svr^t; cos /3 cot (3 (m?^ — m?^) 2 87r2?; cos (3 cot /3 (m?^ — 



""^^ 87r2t; cos P cot /3 (m?^ — m?^) 2 87r2?; cos /3 cot /3 (m?^ — m? 

3m^AA2A^^sin0b (^KV m|) 3m^AA sin log(mym-y 



87r2t) cos /3 cot /5 (m?^ — 771^2 87r2ti cos cot /5 (m?^ — m?^) 
_ _3m^AVAA^^sin06^(my my 



87r2t;2 cos^ P 


(m? • 
^ 02 


- m? )2 


3mlX^xAl^ g{ml, 


m? ) 
62^ 


3mlX' 


Sn'^v cos P cot /3 (m?^ — 


m~ )2 

bi' 


87r2? 


3m^A2xAbA^2 sin 0^ 




' <) 


87r2t; cos P cot /5 


(m? - 

^ 02 


- m? )2 

bi' 


Sm^X'^xAbAi^ sin 0;, 






87r2t; cos2 P cot /? 


(m? - 

^ f)2 


- mi Y 



2 



— m'i \2 8^2,, noQ/^ ■' ^ ^1' ''2 
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where 



Aj^ = At cos (pt + Ax cot (3 , 

Aj^ = At + \x cot P cos 0t , 

A^^ = At + Ax tan /3 cos 06 , 

Ai = y4f,cos</)t + Axtan/3 , 



(13) 



and g{ml,my) is another scale-independent function that is defined as 




(14) 



Among these Mij, those that are responsible for the CP mixing between scalar and 
pseudoscalar Higgs fields are M13, M23, M15, M25, M34, and M45 in the (5*1, S2, Pi, S3, P2) 
basis. They are clearly zero at the tree level. Thus, there is no CP mixing between 
scalar and pseudoscalar neutral Higgs bosons at the tree level. They become non-zero 
as the radiative corrections Mj^ and M^^ are taken into account. Moreover, they receive 
complex phases 0* and (ph from M*^ and M^^, respectively, for the CP mixing. Notice 
that M*3, M*3, M*5, M*5, M*4, and M*5 depend on 0^ while M\^, M^^, M\^, M^^, M^,, 
and M45 depend on 0^. Therefore, the magnitude of the CP mixing between scalar and 
pseudoscalar neutral Higgs bosons is directly dependent on 0^ and 0f,. The CP mixing 
would be maximal when sin0t = sin0b = 1. 

The five eigenvalues of M, denoted as m^, {i = 1-5), define the five physical neutral 
Higgs bosons hi {i = 1-5) as the mass eigenstates and their squared masses. Unless 0^ or 
0b vanishes, hi {i = 1-5) would not have definite CP parities. They are in general given 
as the mixtures of Si, S2, S3, Pi, and P2. We sort these five neutral Higgs bosons in 
the increasing order of their masses such that m^^ is the smallest eigenvalue and hi is 
the lightest neutral Higgs boson. If 0t = 0f, = 0, these five neutral Higgs bosons may be 
classified into three scalar and two pseudoscalar Higgs bosons. 



In the standard model (SM), it is well known that particles acquire their masses through 
the interactions with the Higgs field. This implies the existence of the yet-undiscovered 
SM Higgs boson. The lower bound on the SM Higgs boson mass of about 114.5 GeV is 
determined by means of the negative result for the Higgs searches at LEP2 experiments. 
The SM does not predict its mass. But the SM provides the decay modes as well as 
production rates for each possible mass range. 

The decay signature of the SM Higgs boson depends on its mass. If the Higgs boson 
mass is in the range of 80-140 GeV, the CMS [28] and ATLAS [29] collaborations expect 



III. HIGGS DECAYS 
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that its decay into a pair of photons would be the most interesting channeL In this range, 
other important channels are the decays into pairs of bottom quarks, charmed quarks, 
and tau leptons, as well as a gluon pair. In the mass range between 114.5 GeV and 2 
(=180 GeV), its decays into pairs of weak gauge bosons, where one of the gauge boson 
in the pair is virtual, become dominant besides the decay channel into a bottom quark 
pair. For the SM Higgs boson mass between 2mz and 2mt (= 350 GeV), it would almost 
exclusively decay into the weak gauge boson pairs (IWl^, ZZ). Its mass being larger than 
2mt, the SM Higgs boson would decay mainly into pairs of weak gauge bosons or top 
quarks. 

From the extrapolations of the relevant couplings through the Higgs decay modes in 
our model, the experimental constraint on the mass of the lightest neutral Higgs boson 
in our model might be derived. Therefore, the calculation of the Higgs decay modes is 
essential for the Higgs searches. The total decay width of the jth neutral Higgs boson in 
our model may be assumed as [30] 



T{hj) = T{hj hh) + V{hj rr) + T{hj fifi) + T{hj cc) + T{hj ss) 

+ T{hj gg) + T{hj WW) + T{hj ZZ) + T{hj tt) , (15) 



where the notations for each particle can be understood without difficulty. 

Note that the Higgs bosons might decay into a supersymmetric particle pair if their 
masses are small enough. For a comprehensive analysis, the Higgs decays into sfermions, 
neutralinos, and chargions also must be considered in explicit CP violation scenario on 
the present model. The Higgs decays into a superparticle pair do not play an important 
role for relatively light Higgs bosons [30]. Also, it has been noted that the decay modes 
into supersymmetric particles acquire significant branching ratios and can be dominant 
ones if they are possible [30]. In our case, we set the parameter values such that a common 
SUSY breaking scale is taken to be 1 TeV. Thus, the squark masses are relatively heavy 
in our case, and therefore we speculate that the Higgs decays into a squark pair would 
contribute weakly to the total Higgs decay. We note that our speculation is based on the 
size of the coupling coefficients between Higgs bosons and a pair of squarks, which we 
explicitly derive. Using them, we calculate the partial decay width of a Higgs boson into 
a pair of gluons via squarks. It remains as a future study to calculate the Higgs decays 
into a pair of supersymmetric particles in order to examine our speculations on its relative 
dominance in the total Higgs decay. 

The partial decay width of a neutral Higgs boson hj into a pair of down or up quarks 
is given as 



T{hj — > dd) 



Cfglmlmh^ 
327rm^ ^ 



1 - 





uu 



Cfglmlnih^ 
327rm^ ^ 



1 - 




\^ hjuu) 



( 



1 




) 



+ K.nJ' , (16) 



where 




hjdd 



COS/3 



sin P 
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hiuu 



cot (30. 



3j ' 



(17) 



where Oij are elements of the orthogonal transformation matrix which diagonalizes the 
mass matrix for the five neutral Higgs bosons. The color factor is C/- = 3 for quarks and 
Cf = 1 for leptons. The partial decay width of a neutral Higgs boson hj into a pair of 
charged leptons is the same as T{hj —>■ dd). 

The partial decay width of hj into a pair of weak gauge bosons may be obtained from 
that of the SM Higgs boson, through a relation given by 



Ti^h, VV) = Glyy Tsuih, - VV) 



(18) 



where rsM(^/ VV) is the decay width of the SM Higgs boson into a pair of gauge 
bosons, one of produced gauge bosons being virtual. 



Ghjvv = cosfSOij + sin f302j ■ 
The partial decay width of hj into a gluon pair is given as 



^{hj gg) 



(|A^f +|A^|^) 



(19) 



(20) 



where asijnz) is the coupling coefficient of the strong interactions, evaluated at the elec- 
troweak scale, and and are respectively the scalar and pseudoscalar gluon ampli- 
tudes, which are given as 



A' 
A^ 



E 

q=t,b 

E 

q=t,b 



^hjqq-^qi'^q, 



E ^hjqaqa^_2 ^^(^QJ 



a=l,2 



2m? 

la 



(21) 



where Ghjqaqa is the coupling coefficient of the neutral Higgs bosons to a squark pair, 
Ag{Tq) and (r^) are respectively the scalar and pseudoscalar form factors due to the 
quark, and A^ (r^) is the scalar form factor due to the squark, given as 



A'qir,) 



rjl + (l-r,)/(r,)] , 
with the scaled variables defined as 



A^r. 



^{rq) 



Tn 



Am? 

£6 

2 



Arnj- 
%_ 

2 

mi 



r,-[r,-/(r,-) - 1] , (22) 



(23) 



and the function / defined as 

arcsin^(l/A/r) 

+ 



f{r) 



log 



r 



1% 



r>l, 
r < 1 . 



(24) 



Notice that there is squark contributions to the partial decay width of hj due to the 
couplings of the neutral Higgs boson to a pair of squarks of the third generation. Let us 
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describe the couplings of the neutral Higgs bosons to a squark pair in more detail. They 
are given as 

G.,M™ = M..O.,f/rt/L, (25) 

where is the unitary matrix that transforms the weak eigenstates of the squark q 
{q = t, b) to the mass eigenstates, r'**^'^ is the coupling of hj to a pair of the weak 
eigenstates of the squark q, and the subscript indices run l,m = 1,2 and r, s stand for 
L, R. We may parameterize as 

cos^g -sin^g-e-^'^^ \ 
1^ sin^g-e^^^ cos^^- J ' ^ ^ 

where the mixing angles 6q and the complex phases (pg vary between — 7r/2 and 7r/2, and 
thus both cos dq and cos (pg are non-negative. We obtain the explicit expressions for T^'^^ 
(g = t, b) as follows: 



Xx \ 
\x ) ' 



h-stt 





f ^ sin 




mt 




f 2 sin 




mt 


1 



A cos P \ 



vsin/3 \ Acos/3 
n^f 

I -i(Ae+^<^*cot/3 + Ax) 



j.h,it ^ ^( 2(Ae-***cot/3 + Ax) 



r 



-ph2bb 
-ph^bb 



rrib 


cos 








cos 






1 



iAcot/? y ' 

-yl;,e+*'^' 2mb 
? Ax \ 
1^ Ax j ' 
/ A sin /3 \ 



wcos/5 V Asin/3 



ph466 _ !!!^ f z(Ae-*<^nan/3 + Ax) 



^2 



-z(Afee+*'^''tan/3 + Ax) 



r.... ^ !^ zAtan/3 _ ^^7) 



-zAtan/5 



The above formula for the coupling coefficient between a Higgs boson and a pair of squarks 
may be used for the calculation of the partial decay width of the Higgs boson into a pair 
of stop or sbottom quarks. 



IV. NUMERICAL ANALYSIS 

For our numerical analysis, we fix some parameter values. The renormalization scale is 
taken as A = 300 GeV. The quark masses are taken as mt = 175 GeV for top quark 
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and rrii, = 4 GeV for bottom quark. The strong coupling constant is evaluated at the 
electroweak scale as a{mz) = 0.1187, and the weak mixing angle is set as sin^ 9w = 0.23. 
Then, there are a few free parameters in our model, the MNMSSM with explicit CP 
violation at the one-loop level. The CP symmetry in our model is explicitly violated at 
the one-loop level by the presence of (pt and (pb, arising from the squark contributions of 
the third generation in the effective potential. Other relevant free parameters besides (pt 
and (pb are: tan (3, A, Ax, x, ^, itiq, At, and Ab- We take to be At = Ab, and (pt = (pb for 
simplicity. 

We first calculate rrih-^, the mass of the lightest neutral Higgs boson in our model, 
for tuq = At = 1000 GeV, tan/3 = 5, A = 0.1, Ax = x = 100 GeV. The result is 
shown in Fig. 1, where contours of m/^^ is plotted on the (^, 0j)-plane. Please notice that 
ruh^ increases as the tadpole coefficient increases, because the tadpole term prohibits 
a massless neutral Higgs boson. The (^, 0t)-plane in Fig.l has some shaded regions. 
The dotted region is where the spontaneous symmetry breaking does not work, and the 
hatched regions are experimentally excluded by the LEP2 constraint. Only the remaining 
bright region is physically allowed. Therefore, for the chosen parameter values, we find 
that 12 < rrih^ < 52 GeV. 

The hatched regions, imposed by the LEP2 data, may need some explanations. We 
assume that the SM Higgs boson decays exclusively into a pair of bottom quarks. Then, 
we apply the LEP2 data to obtain the corresponding couphng of hi in our model. We 
interpret the experimental constraint upon the SM Higgs boson by LEP2 for the parameter 
region in our model, with the modified Higgs coupling. In this way, the hatched regions 
are established. 

One may note that the allowed mass range between 12 and 52 GeV for hi in our 
model is significantly below 114.5 GeV, which is the LEP2 lower bound on the SM Higgs 
boson. However, it does not imply in a straightforward way that the present model 
is phenomenologically contradicting the experimental constraint of LEP2. In the CP- 
violating case, the mass of the neutral Higgs boson may be released from the LEP2 
constraint. For example, the MSSM in the explicit CP violation scenario with (pt = 7r/2 
allows the existence of a neutral Higgs boson with a mass as small as 30 GeV [31]. This 
implies that discovering a neutral Higgs boson depends not only on its mass but also 
its coupling coefficients and other factors. Since the lightest neutral Higgs boson in our 
model is a mixture of the neutral components of the three Higgs fields, the possibility of 
its discovery would also be affected by how much the Higgs singlet field is mixed in it. 

We next calculate ruh^ {i = 1-5), the masses of the five neutral Higgs boson in our 
model, for tuq = At = 1000 GeV, tan/3 = 5, A = 0.1, = 100 GeV, x = 200 GeV, and 
C, = 100 GeV. The result is shown in Fig. 2, where m/j. {i = 1-5) are plotted as functions 
of < </)f < TT. It can easily be observed that, for the chosen parameter values, all of the 
five neutral Higgs bosons have small masses, below the top mass. Also, the dependence 
of ruh^ {i = 1-5) on the CP phase (pt is quite recognizable in Fig. 2. As the CP phase 
increases from zero to tt, the figure shows that rrih^ increases from 60 GeV to 66 GeV, 
TJih^ from 68 GeV to 69 GeV, m,,^ from 73 GeV to 94 GeV, ruh^ from 79 GeV to 127 GeV, 
and m/jg from 143 GeV to 163 GeV. Thus, the neutral Higgs bosons deviate 10 %, 1.4 %, 
30 %, 60 %, and 14 % in their masses, respectively, for < 0j < vr. The variation of the 
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neutral Higgs boson masses in the MNMSSM against the CP phase may be regarded to 
be very large as compared to the result in the NMSSM [16], where the five neutral Higgs 
boson masses vary below 5 % against the CP phase. 

We continue our numerical analysis to calculate the branching ratios of each of the five 
neutral Higgs bosons. For the same parameter values as in Fig. 2, we obtain a number of 
important branching ratios. The results are shown in Figs. 3a-3e, where they are plotted 
as functions of < 0t < vr. Note that for the whole range of the CP phase the decay 
channel into a pair of bottom quarks is most dominant for all five neutral Higgs bosons. 
The branching ratio of BR{hi — > bb) is invariably more than 90 % for all of hi {i = 1-5). 

For the branching ratios of h^, the heaviest neutral Higgs boson in our model, are 
somewhat different from the other four neutral Higgs bosons. Since its mass is calculated 
to be between 143 GeV and 163 GeV for < 0( < tt, its decay into a pair of weak gauge 
boson might be important, as Fig. 3e indicates. It is noticeable that BR{h^ WW) is 
comparable to BR{h^ —>■ bb) for (pt ~ vr. The fluctuations of the branching ratios of 
against the CP phase are also recognizable in Fig. 3e. For decay, BR{h^ t^t~) and 
BR{h^ ZZ) remains almost stable, whereas BR{h^ cc) and BR{h^ gg) decrease 
by an order of magnitude, as (pt varies from zero to tt. 

The patterns of the branching ratios of hi, h2, h^, and /i4, are nearly similar to each 
other. The next dominant decay channel is hi t^t~ for i = 1-4. The branching ratios 
do not wildly fluctuate against the CP phase, and the relative size between the branching 
ratios are BRQii bb) > BRQii r+r") > BRQii — > cc) > BR{hi gg), for i = 1-4. 

We would like to remark the interesting behavior of h^ at the branching ratio for the 
gluon channel BRiJi^ gg). As shown in Fig. 3d, we have BR{h4 —>■ gg) = 0.9852 x 10~^ 
for (pt ~ 0.9331. We note that the most part of this value comes from the scalar gluon 
amplitude rather than the pseudoscalar gluon amplitude for (pt = 0.9331. This is 
because the top and bottom quark contributions cancel out accidently each other in the 
pseudoscalar amplitude A^ for h4^ in Eq. (20), That is, Gf^^^Af (r^) and Gf^j^f^Af (rt) in A^ 
have accidentally almost the same value but with opposite sign, for (pt = 0.9331. Thus, 
the scalar gluon amplitude contributes dominantly to the BR{hi gg) in our analysis. 

Finally, we calculate T{hi), the total decay widths of hi {i = 1-5), for the same param- 
eter values as in Fig. 2. The result is shown in Fig. 4, where T{hi) {i = 1-5) are plotted 
as functions of < 0( < tt. Their dependence on the CP phase is quite large. While 
T{hi) and r(/i2) decrease as (pt increases, the total decay widths of the other neutral Higgs 
bosons increase with increasing (pt. In particular, it is worthwhile noticing the increase 
of r(/i5). This behavior is understandable because the decay channel into a pair of weak 
gauge bosons is only allowed for h^. 

V. CONCLUSIONS 

We study the Higgs sector of the MNMSSM with tadpole terms. We see that the model 
can accommodate explicit CP violation at the one-loop level. As the squarks of the third 
generation with non-degenerate masses give rise to the radiative corrections, the explicit 
CP violation can be generated by complex phases (pt and (pf, which appear in their mass 



14 



matrices. At the tree level, neither explicit nor spontaneous CP violation would be viable 
in the Higgs sector of the model. Unlike our model, the NMSSM may possess complex 
phases at the tree level. 

In the presence of the complex phases that trigger the explicit CP violation at the 
one-loop level, we calculate the masses, the branching ratios for dominant decay channels, 
and the decay widths of the five neutral Higgs bosons in our model. The masses show 
fluctuations up to 60 % as the CP phase varies, for given parameter values. Whereas 
the heaviest neutral Higgs boson in our model can be as heavy as 160 GeV, the rest of 
them is relatively light. In particular, the mass of the lightest neutral Higgs boson in our 
model is predicted to be as small as 12 GeV, for reasonable ranges of relevant parameter 
values. Although the predicted mass of the lightest neutral Higgs boson is well below 
experimental constraint of LEP2, it might have been escaped from experiments since its 
discovery depends on its coupling coefficients to other particles as well as on its mass. 

The coupling coefficients of the neutral Higgs bosons to fermion pairs, weak gauge 
boson pairs, gluon pairs, and squark pairs have been explicitly calculated. The results 
of calculation, which have not been derived before, may be used helpfully for the search 
of Higgs boson in similar supersymmetric models in the future experiments. In terms 
of these coupling coefficients, the branching ratios of the five neutral Higgs bosons for 
dominant decay channels are calculated, where the CP violation effect is included. We 
find that the heaviest neutral Higgs boson exhibits fluctuating branching ratios against 
the variation of the CP phase. For the rest four Higgs bosons, the decay into a pair of 
bottom quarks is most dominant, and the branching ratios are rather stable against the 
variation of the CP phase. Finally, we calculate the total decay widths of the five neutral 
Higgs bosons, which are significantly dependent on the CP phase. 

If CP is violated at the one-loop level in the MNMSSM with tadpole terms, the result 
of our calculations, which has been done for a representative set or range of parameter 
values, suggests that the whole parameter space should be investigated in order to draw 
more interesting observations within the context of CP violation. 
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FIGURE CAPTIONS 



Fig. 1: Contours of the lightest neutral Higgs boson mass on (/)t)-plane, for niQ = At = 
1000 GeV, tan/5 = 5, A = 0.1, and Ax = x = 100 GeV. The bright region is physically 
allowed whereas the shaded regions are excluded. The spontaneous symmetry breaking 
does not work in the dotted region, and the experimental constraint by the LEP2 data 
excludes the hatched regions. 

Fig. 2: The masses of the five neutral Higgs bosons as functions of where the pa- 
rameters are set as mq = At = 1000 GeV, tan/5 = 5, A = 0.1, Ax = x/2 = ^ = 100 
GeV. 

Fig. 3(a): The branching ratios of hi as functions of (pt, for the same parameter values as 
Fig. 2. 

Fig. 3(b): The branching ratios of /12 as functions of 0^, for the same parameter values 
as Fig. 2. 

Fig. 3(c): The branching ratios of /13 as functions of (pt, for the same parameter values as 
Fig. 2. 

Fig. 3(d): The branching ratios of as functions of (pt, for the same parameter values 
as Fig. 2. 

Fig. 3(e): The branching ratios of /15 as functions of (pt, for the same parameter values as 
Fig. 2. 

Fig. 4: The total decay widths of the five neutral Higgs bosons as functions of (pt, for the 
same parameter values as Fig. 2. 
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Figure 1: Contours of the lightest neutral Higgs boson mass on (,^, </)t)-plane, for mq = 
At = 1000 GeV, tan/5 = 5, A = 0.1, and Ax = x = 100 GeV. The bright region is 
physically allowed whereas the shaded regions are excluded. The spontaneous symmetry 
breaking does not work in the dotted region, and the experimental constraint by the LEP2 
data excludes the hatched regions. 
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Figure 2: The masses of the five neutral Higgs bosons as functions of 0t, where the 
parameters are set as mq = At = 1000 GeV, tan/5 = 5, A = 0.1, A\ = x/2 = C, = 100 
GeV. 
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Figure 3a: The branching ratios of hi as functions of (pt, for the same parameters as Fij 
2. 
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Figure 3b: The branching ratios of /i2 as functions of (pt, for the same parameters as Fij 
2. 
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Figure 3c: The branching ratios of /13 as functions of (pt, for the same parameters as Fij 
2. 
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Figure 3d: The branching ratios of /i4 as functions of (pt, for the same parameters as Fij 
2. 
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Figure 4: The total decay widths of the five neutral Higgs bosons as functions of (pt, for 
the same parameters as Fig. 2. 
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